Cavity-assisted energy relaxation for quantum many-body simulations 
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We propose an energy relaxation mechanism whereby strongly correlated spin systems decay into 
their ground states. The relaxation is driven by cavity quantum electrodynamics interaction and the 
decay of cavity photons. It is shown that by applying broadband driving fields, strongly correlated 
systems can be cooled regardless of the specific details of their energy level profile. The scheme 
would also have significant implications in other contexts, such as adiabatic quantum computation 
and steady-state entanglement in dissipative systems. 



Introduction. — Advances in laser cooling and trapping 
technologies have opened up intriguing possibilities of us- 
ing atomic systems to simulate quantum many-body phe- 
nomena in strongly correlated regimes Concerning 
such quantum simulation, diverse studies have been car- 
ried out both experimentally and theoretically, address- 
ing various fundamental many-body problems Q ■ From 
the theoretical perspective, the properties of the ground 
states are among the main concerns regarding strongly 
correlated systems. On the other hand, studies on the 
experimental realizations have rather placed the main 
emphasis on how to mimic particular many-body model 
Hamiltonians, while it is often overlooked that being able 
to realize a particular Hamiltonian does not mean that 
one can prepare its energy states. 

The most natural way to bridge such a logical gap 
would be to introduce energy relaxation in simulated 
strongly correlated systems. First of all, it is impor- 
tant to understand that in strongly correlated regimes 
no conventional cooling scheme provides a proper relax- 
ation channel 0, Q]. The main reason is because they 
were devised to cool weakly correlated systems, keeping 
it in mind that cooling individual atoms results in cool- 
ing the whole system. This does not hold in strongly 
correlated regimes. Moreover, while important classes 
of many-body systems emerge when the atomic internal 
degrees of freedom (DOFs), i.e., spins, are strongly cor- 
related, cooling of such systems is yet to be considered. 
In weakly correlated regimes, optical pumping could be 
regarded as the spin counterpart of the conventional cool- 
ing of atomic motions. However, the problem again be- 
comes elusive when the spins are strongly correlated. 

Up until now, the adiabatic method is the only avail- 
able way to prepare the ground state of a given strongly 
correlated system The underlying idea is to prepare 
a known ground state in a weakly correlated regime and 
transfer it to that in a desired strongly correlated regime 
by adiabatically changing the system parameter. This 
method is feasible when the ground-state energy gap is 
finite throughout the whole process. However, this condi- 
tion is difficult to achieve because such a transition (from 
a weakly correlated to a strongly correlated regime) is 
usually a quantum phase transition, at the critical point 



of which the energy gap vanishes in the thermodynamic 
limit [5(. Moreover, this method requires some a pri- 
ori knowledge about the system, while quantum simula- 
tion would be more useful when one has little informa- 
tion about it. As adiabatic dynamics is a slow unitary 
evolution, the coherence time of the system also limits 
its performance both during and after the preparation, 
whereas energy relaxation would continuously stabilize 
the ground state competing with decoherence. 

In this paper, we introduce a cavity-assisted energy 
relaxation mechanism that directly cools strongly corre- 
lated atomic systems. As a concrete example, we discuss 
cooling of strongly correlated spin systems [a, LD, [8( • Al- 
though some schemes have been proposed to prepare par- 
ticular ground states 0, [l(| , there has hitherto been no 
genuine idea of energy relaxation applicable to general 
many-body Hamiltonians whose ground states are un- 
known. Besides the application in quantum simulations, 
the importance of energy relaxation can also be found in 
many other contexts. For example, it has a potential to 
provide a more efficient route to reach the ground state 
for adiabatic quantum computation 11]. Also, recall- 
ing that strongly correlated spin systems have entangled 
ground states, it would provide a general framework for 
establishing steady-state entanglement, discussed so far 
from different perspectives [10j. Our scheme is particu- 
larly timely amid huge current interests in placing many- 
body atomic systems inside cavities 12| with some real 
experiments [13| . 

Physical model. — We consider strongly correlated spin 
systems as can be simulated in various atomic systems, 
such as optical lattices @, ion traps 0], and coupled 
cavities 8] . The common feature of those quantum sim- 
ulations is that atomic hyperfine ground levels represent 
spins, while some external DOF (atomic momentum, cav- 
ity modes, etc.) mediates their coupling. In so doing, the 
external DOF is effectively decoupled from the dynam- 
ics. Our concern is how to manipulate the spin DOF to 
lower its respective energy. 

For simplicity, we present our idea with s = h spin 
systems. Two ground levels of an atom represent spin- 
down and up states, as shown in Fig. [TJ To begin with, 
let us write the Hamiltonian of the quantum simulator as 
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H o = J2^=o E n where is the /xth eigenen- 

ergy with < E v for « < z/, is the /zth eigen- 

state, and TV is the number of spins. We assume that 
the ground state \*&q) is nondegenerate. Our aim is to 
figure out a process that allows cooling (l^) — > \^f v ) 
for fi > v) but forbids heating (« < v). First, we 
note that in the aforementioned simulations, spontaneous 
emission is not feasible as a relaxation channel for three 
main reasons: (i) As the spontaneous emission rate is 
larger than the characteristic energy scale of the simu- 
lator, the energy levels can not be resolved due to the 
line broadening, (ii) The spontaneous emission destroys 
the atomic coherence, hence the coherences within the 
many-body states, (iii) As atomic excited levels are out 
of the spin subspace in consideration, the energy of the 
simulator is not clearly defined when the atom is ex- 
cited. For these reasons, we need to introduce some 
ancillary DOF so that energy is first transferred coher- 
ently to it and then relaxed from it. Let us denote by 
l^^n) the state of the system along with the ancillary 
mode in the nth excited state. Our strategy is to find 
such a transition as |\& «-> \^„,n+ 1} — > |^„,n) for 
u > v. Note that the irreversible transition in the inter- 
nal DOF is attained by destroying the coherence of the 
external DOF. It is worth comparing this with the side- 
band cooling that works in an opposite way. Denoting 
by g) and |e) the atomic ground and the excited states, 
respectively, the sideband cooling can be represented as 
\g,n) <-> |e,n- 1) -> \g,n- 1). 

It can be seen that the cavity field is best suited for 
the above purpose. Among various possibilities, we con- 
sider a situation where one cavity is coupled to the quan- 
tum simulator, as depicted in Fig. Q] Two transitions 
\ e )j an d lt)j \ e )j °f the jth atom are cou- 
pled, respectively, to the two orthogonally polarized cav- 
ity modes, whose annihilation operators are denoted by 
a\ and a 2 , with coupling rates g±j and g 2 j and detunings 
Ai and A2. The classical field with Rabi frequency £l X j 
{x = 1, 2) is also applied with detuning A x + S X . We take 
a large detuning regime where A^, Ai — A 2 ^>g X j, &xj, S x . 
In this regime, the atomic excitation is suppressed and 
the conventional adiabatic elimination method can be 
used to yield the effective Hamiltonian [3] 

H\ = Hd + H t , 

H d = ~e z - eia\ai - e 2 a\a 2 , (1) 
H t = -{a\T-e i6lt + a\Y + e iS2t + H.c), 

where e z = £j _ Jl^) s * £l = £. l^sjsf, 

= 1 = Ei^v and 1 = 

Y^j 9 \ 2 2j Here, we employed the spin notation with 
s z = f(\V (T|-U)ai)and S + = ( S -)t = |t)a|. 

Working principle. — Let us take two of the energy 
eigenstates j^} and j'S,,) of the spin system. In what fol- 
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FIG. 1: Atomic levels and transitions. Two transitions |j) . <-» 
\e) ■ and \]) ■ <-> \e) ■ are coupled, respectively, to the two 
orthogonally polarized cavity modes. Two Raman transitions 
represented by the subscripts 1 and 2 are embedded. 

lows, we adopt the convention that (W^J O \^ v ) for any 
operator O is denoted by {0)^ u . Our plan is to treat 
the Hamiltonian |T]) as a time-dependent perturbation. 
For this, we impose |25 M „| > (e 2 ) w , (e z ) vv , |(e z ) M „|, and 
similarly for ei, £2, T±, and T±, where E^ = — E v . 

The Hamiltonian (p} comprises the diagonal term Hd 
and the off-diagonal term H t . The former makes a small 
correction to the energy levels, which we take into ac- 
count later. The essential part is the latter that induces 
transition between the energy levels. Let us denote by 
|\E , M ,n 1 ,n 2 ) the state with n\ and n 2 photons in the ai 
and a 2 modes, respectively. If the cavity has initially 
no photon, the transition occurs creating a cavity pho- 
ton. The transition rate is thus determined by the matrix 
elements 1, 0| Hi 0, 0) = {T-) fll ,e^ E ^ +Sl '> t and 
(1V,0,l|ff, |*„,0,0> = (r + )^e l ( E -+ fe ) t , where H T = 
e iH tjj te -zH t j g foe interaction Hamiltonian. This ex- 
pression indicates that the transition occurs dominantly 
when the resonance condition E^ v + <5{i.2} = is satis- 
fied. Furthermore, if ^{1,2} is chosen to be positive, the 
resonance condition is satisfied only when u < i>, which 
means that the transition occurs in favor of decreasing 
the energy of the quantum simulator. <5{i,2} being posi- 
tive implies that the two-photon Raman transition takes 
place in such a way that an atom absorbs a lower-energy 
photon from the classical field and emits a higher-energy 
photon into the cavity mode. The net result is thus 
that the amount of energy <5{i,2} is transferred from the 
quantum simulator to the cavity mode. The subsequent 
decay of the cavity photon then completes one cycle of 
the energy relaxation. The correction by Hd in Hamilto- 
nian (TTJ, albeit small, can be taken into account as fol- 
lows. Hd shifts the energy of state l 1 !^, Tii, n 2 ) as much 
as — (e z + ?ii£i + n 2 e 2 )nn, which modifies the resonance 
condition as E^ v - (e z + £{i, 2 } W + ( e zW + ^{1,2} = 0. 
This term also modifies the energy eigenstates them- 
selves. The infidelity of the state for n\ = n 2 = due to 
this correction is of order \(e z ) fl ^/E^\ 2 . 

We impose the following additional conditions. First, 
we require \E^ V \ ^> k, where k is the cavity decay 
rate. Otherwise, the energy levels can not be resolved 
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FIG. 2: Numerical simulation for a two-spin system, (a) 
Time evolution of the population (top panel) and detection 
rates of the cavity photons (middle and bottom pannels). 
(b) The time evolution of the ground state population in 
the presence of the thermal field with mean photon number 
n = {0, 0.02, 0.04, 0.06, 0.08, 0.1} (from top to bottom). 



due to the line broadening, which results in heating 
of the quantum simulator. Second, for the lower en- 
ergy state, (r±) MM <C K should be satisfied, since oth- 
erwise the cavity photon can be generated without a 
transition, resulting in heating, e.g., through transition 
1*^,1,0) -> 1^,0,0) that is allowed for fj, < v. Third, 
for a similar reason, there should be no thermal photon 
coming into the cavity from the environment. 

As an example, we analyze in Fig. [5] the energy relax- 
ation process of a two s — ^ spin system with Hamilto- 
nian Hq = B(s{ + sf) + Jsi ■ s 2 and B = 2J > 0. In 
this case -E^+i — E^ = B for every fj,. The parameters 
are chosen as A/y^n/lOO = Oy/lO = gtijzj} = ^ K anci 
^{1,2} = 10k = B. H,2j is chosen to be —Clij to optimize 
the process. e z is taken to be zero. The top panel of 
Fig. fJK shows the population of each energy level with 
respect to time and the middle and the bottom panels 
show, respectively, the detection rates of the photons in 
mode a\ and a 2 averaged over the correspon ding time 
bins by using the quantum trajectory method [15]. The 
population, initialized in the highest energy level to show 
the cascade-like transition structure, eventually trapped 
in the ground level, during which the detection eventually 
ceases, indicating the completion of the relaxation. We 
show in Fig. [Jb the ground state population for various 
mean thermal photon number n. As discussed before, 
the asymptotic fidelity decreases as n increases. 

Note that this simple model already illustrates a novel 
application of the present mechanism. The ground state, 
which is a maximally entangled state, is prepared as a 
steady state without any initialization of the system or 
time control of the fields. For a small number of atoms, 
one can apply several resonant fields in the same man- 
ner to prepare the entangled ground state. The effective 
spontaneous emission rate ~ Xj -7|f2j/A| 2 , with 7 the 
intrinsic spontaneous emission rate, can be made arbi- 
trarily negligible compared to the Raman transition fre- 
quency ~ ■ \gjflj/A\ by increasing A. 

General scheme. — The above analysis leads us to the 
reasoning that if the classical fields are broadband with 



Si t 2 > 0, many-body systems would tend to cool even if 
they have many energy levels with different energy gaps, 
hence many different resonance conditions. Let us first 
simplify the situation by considering spectrally incoher- 
ent broadband fields. We denote the spectral density 
of the broadband field by 1(5), which is normalized as 
J I(6)dS = 1 (the subscripts 1 and 2 are omitted for 
brevity). The spectral incoherence can be reflected by 
replacing e m with J d8f(8)e lSt and J dSdd' f(5)f*{5') ex- 
plicitly with J ddd6'I(6)d(6 — 6'), where X>(x) is the Dirac 
delta function. The standard treatment of the time- 
dependent perturbation theory then yields the Fermi- 
golden-rule-like transition rate 27rni |(r_) Mt/ | 2 I(E VI1 ) f° r 
transition n\ — 1, n 2 ) — > n 2 ) (similarly for 

n 2)- Hd is invariant under the above change, making the 
same contribution. 

Now that the transition is understood to be Marko- 
vian, it can be represented by a Markov chain, where 
each node represents an energy level (including the cav- 
ity state). For each pair of the energy levels, the tran- 
sition is established if the broadband field contains the 
spectral component 8 corresponding to the energy gap. 
The cavity decay |\t ^ni,^) — ► \^^,rii — 1,712) is also 
included with rate n\K (similarly for n 2 ). 

The detailed relaxation process depends on the spe- 
cific many-body systems at hand. In what follows, we 
discuss the relaxation process of the isotropic Heisen- 
berg spin chain with Hamiltonian Hq = J X^Li* $j ' > 
where J > and N is even for \^q) to be nondegener- 
ate. Those arguments below, however, mostly hold when 
Hq commutes with the total spin operator S z = J^j s j; 
as is the case for any Heisenberg spin system without 
XY anisotropy. Most of the spin systems simulable by 
earlier schemes fall into this class. Several remarks are 
in order, (i) All the conditions mentioned above are in- 
herited. However, we strictly impose them only between 
^o) and every with /i ^ 0, since this is important 
for the ground state to remain stable. Even if the con- 
ditions are not met for some higher levels, the overall 
effect still pushes the population downwards and once it 
reaches the ground state, any further transition is sup- 
pressed, (ii) we lift most degeneracies by applying a mag- 
netic field BS Z , which can be done by an additional field 
along with the e z term. As [Hq,S z ], the eigenstates of 
Hq are now the simultaneous eigenstates of S z . |Wo) re- 
mains unchanged, (iii) we break the spatial coherence 
by applying spatially incoherent fields, or alternatively 
by breaking the symmetry, e.g., by tiling, focusing, or 
displacing the fields and applying those complementary 
field configurations in turn. In what follows, we simply 
assume |(r ± )^| 2 = ^ |^( S ±U| 2 . 

Efficiency. — The dependency of the relaxation time on 
N is an important issue. If it grows exponentially with 
iV", the relaxation scheme would be much less useful. Al- 
though it is daunting to obtain a general relationship, 
one can see that it grows at most polynomially with TV 
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FIG. 3: Numerical simulation for Heisenberg spin chains, (a) 
The asymptotic population in the ground state, (b) The time 
evolution of the ground state population for gj = 40k. 



mixed state for gj = 40k. Note that as explained above, 
the spectral incoherence, introduced just for ease of anal- 
ysis, makes the process inefficient, which is why such a 
strong atom-cavity coupling is required. 
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for gapped systems as explained below (note that gapped 
systems are usually of special interest). The relaxation 
time can be estimated as the number of energy levels to 
pass through (in the Markov chain) to reach the ground 
level, say iVg, multiplied by the characteristic transition 
time ~ 1/2tt\(T±)^ u \ 2 I(S) between two levels. Although 
the coefficient \(r±) fi v\ 2 is bounded by Eio, it can be kept 
finite for gapped systems. Besides this, we end up with 
a trade-off problem: as the fields get broader-band, Ns 
decreases while 1/1(5) increases. By noticing that the 
eigenstates of Ho are also the eigenstates of S z (ranging 
from —N/2 to N/2) and a transition changes S z by one, 
we see that if the broadband covers the entire energy 
spectrum, Ns is at most 2N (N for spin flips and N for 
cavity decays), while 1/1(5) behaves as N . 

Spontaneous emission. — Although the atomic excita- 
tion is highly suppressed, the spontaneous emission is the 
prominent source of heating. The spontaneous emission 
changes the state as >F) — > s^ or s^sj |\E'), which can 
be reflected in the Markov chain by adding transitions 
|*„> «-> I*,} with rates 3 (l I *)\^ / A\ 2 [\(s+)^\ 2 + 
\( s J)fiv\ 2 + 2|(s|) /J „| 2 ] ! where the branching ratio is as- 
sumed to be 50 : 50. Note that when spectrally incoher- 
ent broadband fields are used, the transition rate is of 
the same order ~ Ifij/Aj 2 as the spontaneous emission. 
In such a case, strong atom-cavity coupling gj ^> n is re- 
quired to obtain a high fidelity to the ground state. See 
the numerical simulation section below. On the other 
hand, when spectrally coherent broadband fields, i.e., 
pulses, are used, the spontaneous emission again can be 
made negligible, in principle, even for weak atom-cavity 
coupling by increasing A, as in the previous case of Fig. [5] 
If the spin system is realized with trapped atoms, the 
spontaneous emission may also excite the atomic momen- 
tum. This could be overcome by sympathetic cooling Q • 

Numerical simulation. — In Fig. we plot the asymp- 
totic ground state population with respect to gj for 
N = {4,6,8}. We take A/gj = g /n and k = 7 = B/10 
with B — Eio/2. Oj is chosen so that | (r+ ) 10 1 = K - 
is constant for 0.5B < 5 — (ei)oo < 3.51? and zero other- 
wise. The spontaneous emission is included as explained 
below. In Fig. [5Jd, we plot the time evolution of the pop- 
ulation in the ground state starting from the maximally 
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